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Abstract 

We present a universal scheme of pulsed operations for the IBM oscillator-stabilized flux qubit com- 
prising the controlled-cTz (cphase) gate, single-qubit preparations and measurements. Based on numerical 
simulations, we argue that the error rates for these operations can be as low as about .5% and that noise 
is highly biased, with phase errors being stronger than all other types of errors by a factor of nearly 10^. 
In contrast, the design of a controlled- ctx (cnot) gate for this system with an error rate of less than 
about 1.2% seems extremely challenging. We propose a special encoding which exploits the noise bias 
allowing us to implement a logical cnot gate where phase errors and all other types of errors have nearly 
balanced rates of about .4%. Our results illustrate how the design of an encoding scheme can be adjusted 
and optimized according to the available physical operations and the particular noise characteristics of 
experimental devices. 

After years of painstaking labor, superconducting qubits [1 are taking shape as viable elements for the 
construction of a scalable quantum computer. Since the initial demonstrations of coherent quantum dynamics 
in superconducting qubits ^[4l[3l|5], it has been recognized that these systems have great potential versatility 
m [3 E], so that one can genuinely envision a quantum-computing integrated circuit emerging from this 
research. However, no clear way forward has been announced, owing largely to one undeniable feature of 
large-scale quantum computation: it will require a very high degree of fidelity in the execution of quantum 
operations, much higher than has been reported in any present experiments. 

How high a fidelity, or how low an error rate, will be needed? On the basis of fundamental early theoretical 
work [HITnilll], lip service is frequently paid to a necessary universal set of operations containing the two- 
qubit controlled-cTx (cnot) gate, and a necessary "threshold" error rate in the 10~^ — 10~^ range. Some 
recent modeling for superconducting qubits [T2j |T3] suggests that such noise levels could conceivably be 
reached in the lab; however, in current experimental practice the ability even to reliably detect such small 
error rates, let alone to achieve them, is in fact very questionable. 

In this paper, we will consider the possibility of constructing a universal set of operations for the IBM 
"Koch qubit" [TT, "8^ Q Although our set of operations will not contain the cnot gate, we will propose an 
encoding scheme for implementing logical cnot gates which can then be used for fault-tolerant quantum 
computation. With the combination of our encoding scheme and other improvements in the theory of 
quantum fault-tolerance [l3|l6l|l3, error rates for our elementary operations in the .1%— .5% range are 
expected to be tolerable for practical quantum computation. We estimate using numerical simulations that 
error rates in this range are possible for the Koch qubit. Though these estimates are far from the error 
rates that have been currently measured in experiments, we hope our results will motivate and stimulate 
the research in superconducting qubits in the IBM lab and elsewhere. More generally and even beyond 
superconducting qubits, our encoding scheme illustrates how techniques of quantum error correction and 
fault tolerance can be tailored to the available physical operations and the particular noise characteristics of 
experimental devices. 

*IBM Watson Research Center, Yorktown Heights, NY 10598. 
^ These authors contributed equally to this work. 

■1- Institute for Quantum Information, California Institute of Technology, CA 91125. 

-•^ After Roger Hilsen Koch (1950-2007), the leader of the experimental superconducting-qubit effort at IBM Research until 
his sudden death on August 4, 2007. 
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Figure 1: On the left, a cnot gate is applied between two qubits. In the middle, the cnot gate is 
simulated by using one cphase and two H gates. On the right, another simulation where the two H gates 
are implemented by "teleportation" using ancilla qubits and measurements. Vertical lines with dots on 
both ends denote CPHASE gates. A triangle pointing to the right denotes the preparation of a qubit in the 
state |+). A triangle pointing to the left denotes a measurement of a^. Conditioned on the measurement 
outcomes, correction operators may be necessary as shown. 



Our proposal involves a synthesis of recent experimental and theoretical developments. On the experimen- 
tal side [TT, 'F, pulsed operations for the Koch qubit were discussed in [18 ; the set of operations considered 
in ^ included the controlled-dz (cphase) gate, the Hadamard (H) gate, single-qubit "diagonal" rotations 
of the form exp(z^crz), the preparation of a qubit in the state |+) = -^(lO) + |1)), and the measurement of 
cTx (where dx, cry, and ciz are the Pauli spin operators). For these operations, \1S\ numerically estimated 
the error rates by considering all physical sources of noise that we presently know how to model for this 
system — 1/f flux noise, instrumental jitter in pulse timing and amplitude, and Johnson noise from resis- 
tances in the circuit. For the cphase gate it estimated an error rate of about .45%, for the H gate about 
.4%, for diagonal rotations about .001% for any ^, for the preparation of the state |+) about .3%, and for 
the measurement of about .2%0 

In contrast, no direct implementation of a cnot gate with error rate close to even 5% could be devised for 
the Koch qubit. Alternatively, a cnot gate can be implemented indirectly by composing other elementary 
operations. Fig. [l] shows two possible methods; the first indirect implementation uses one cphase and two 
H gates, while the second one uses three cphase gates, two qubits prepared in the state |+), and two 
measurements of a^. Therefore, we estimate that the first indirect implementation would have an error rate 
of about 1.25%, and the second one about 2.3%. 

Based on the modeling in [18] , it was further observed that noise in the cphase gate acted predominantly 
as dephasing; i.e., it affected the relative phases in the wavefunction giving rise to "phase" errors which 
can be expressed as diagonal matrices in the computational basis. All other errors, including errors due to 
relaxation and also "leakage" errors associated with transitions to states outside the computational space, 
were observed to be about an order of magnitude weaker. The same observation also applies to single-qubit 
diagonal rotations but, on the other hand, no similar bias was observed for the noise in the H gate. 

On the theoretical side, the findings of [18] motivated us to revisit a longstanding problem in the theory of 
fault-tolerant quantum computation: how to formulate an encoding scheme that exploits large asymmetries 
in the structure of noise in elementary operations. Here and in our companion paper [19 , we present such a 
scheme based on encoding the noisy qubits using a repetition code, and we show how fault-tolerant logical 
CNOT gates can be implemented for this code. The intuition we gained from [18 is that highly biased noise, 
with phase errors being much stronger than all other types of errors, is possible for gates such as the CPHASE, 
but it should not be physically expected for generic operations. Furthermore, an encoding scheme has to 
address the problem that a noise bias can easily be lost as elementary operations are composed together. 

Our solution is to choose a universal set of elementary operations whose implementation induces noise 
that is biased towards dephasing, and where all gates commute with so that the noise bias is maintained. 
The operations we will use are the preparation of the state |+), the controlled-dz (cphase) gate, and 
measurements of observables in the equator of the Bloch sphere, i.e., of the form ex.p{iOaz)cr^ for certain 
angles]^ Noise for the preparation of the state |+) is naturally biased since |+) is an eigenstate of cr^. The 

These numbers are obtained by using the results in [18| and the definition of an error rate in the Supplementary Material. 
In |18| , the entanglement fidelity and not the error rate is used as a measure of the noise strength. 

^Alternatively, we may restrict to measurements of ctx, if we add to our set of elementary operations the preparation of the 
states = exp(i^crz)|+} and \T) = exp(z|^crz)|+}; this is the universal set considered in 19 . For the Koch qubit, it is more 
natural to move all rotations of the form exp(i^crz) before the measurements since, as we will discuss, corrective single-qubit 
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structure of noise for the cphase gate depends on the physical implementation and it can be engineered 
to be biased; here, we propose such a biased-noise implementation for the Koch qubit based on adiabatic 
pulses so that transitions between the computational-basis states are strongly suppressed. Finally, noise in a 
single-qubit measurement can be described by errors preceding the ideal measurement which need not have 
any specific structure since measurement has a classical output. 

These theoretical ideas may have implications for the experiments with the Koch qubit, and solid-state 
qubits in general. They suggest that a profitable focus of experiments could be to design qubits with long 
relaxation time Ti. Provided this is achieved, dephasing, which is the dominant source of noise, can be much 
more effectively suppressed by using our encoding scheme. Furthermore, the implementation of low-noise 
CNOT or H gates is not necessary; it suffices to implement cphase gates with highly biased noise, together 
with single-qubit preparations and measurements. For the Koch qubit, eliminating the need to implement 
the H gate allows us to re-examine the pulse schemes in [18 , and to find a simpler implementation of the 
CPHASE gate for which noise is much more biased than in [18 . Similar simplifications may be possible in other 
types of superconducting qubits. In particular, a standard experimental approach to suppressing dephasing 
noise is to intersperse quantum gates with additional corrective "hardware" operations such as spin echo 
pulses. Our encoding scheme can be seen as a complementary solution where the hardware operations are 
restricted to a much smaller set and dephasing is suppressed at the "software" level by using error correction 
applied on blocks of several noisy qubits. 

We will now discuss the details of our proposal. The central question is how to devise a strategy for 
fault-tolerant quantum computation that effectively exploits a bias in the noise of the cphase gate. To 
develop some intuition, let us assume for the moment that there is no leakage, that independent phase errors 
occur with probability £, and that all other types of errors occur with probability e' <C e. Then a clear 
strategy for the lowest-level coding of quantum data is to use an n-qubit repetition code. One logical qubit 
is encoded in a "block" of n physical qubits, a logical |+) state is represented by all n qubits being in the 
state |+), and a logical |— ) by all n qubits being in the state | — ). Since the logical ciz operator is a dz acting 
on all n qubits in the block, phase errors on more than half of the qubits in the block are necessary for a 
logical error to occur after error correction; taking n odd, and assuming that for each qubit there are t time 
steps where a phase error may occur, the probability of a logical error is approximately 

^L^l^^)ite)'^ . (1) 

On the other hand, the logical <Jx operator is a <Jx acting on any qubit in the block, so that even a single error 
different than a phase error on any qubit and at any time step cannot be corrected; thus, the probability of 
a logical error is approximately 

SL^nte' . (2) 

If the noise "bias" e/e' is large, we can choose some large n and obtain a significant reduction of the logical 
error rate; e.^., if the bias is 10^ and even for t- s = 5%, we find that ~ < 3.5 x 10~^ by setting n = 7. 
This sends an optimistic message, since several schemes are known for effectively implementing fault-tolerant 
quantum computation with error rates of order 10~^ [20 . 

As we will discuss below, our pulsed implementation of a cphase gate can in fact be optimized to lead to 
a noise bias of order 10^. But first, we must examine more closely whether there is actually a gain in using 
the repetition code as our naive argument so far indicates. The main concern is that a large intrinsic noise 
asymmetry in our elementary operations might be spoiled as these operations are composed together. 

For example, consider the problem of implementing a logical CNOT gate between two blocks of encoded 
qubits. The logical CNOT gate for the repetition code can be implemented by bitwise CNOT gates, and each 
CNOT gate can be simulated by using cphase gates as in Fig. [l] But then due to the H gates, phase errors 
that occur during the simulation of each CNOT gate are converted into errors of other types; e.g.^ a <Jz error 
during the implementation of a H gate will be converted to some linear combination of a dz, a dx, and a 
cTy error. So, this construction of a logical CNOT gate destroys the asymmetric structure of the noise and 
nullifies the effectiveness of the repetition code. Avoiding this interconversion of noise is possible, but it 
requires a circuit of greater complexity. At the heart of this circuit construction is the identity operation in 

diagonal rotations are necessary in the implementation of every CPHASE gate. 
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Figure 2: (a) On the left, an identity circuit which consists of preparing a qubit in the state |+), measuring 
(Jz (Jx on this qubit and the input qubit, and finahy measuring on the input qubit. After applying 
the correction operators shown conditioned on the measurement outcomes, the state of the input qubit is 
"teleported" to the output qubit. On the right, the measurement of cj^ cj^ is implemented by using an 
"ancilla" qubit (shown in green) which interacts with the other two "data" qubits (shown in red) via CPHASE 
gates before it is measured, (b) The circuit in (a) where the data qubits are encoded in the n-bit repetition 
code (here n = 3); the ancilla qubit remains unencoded, but the subcircuit enclosed by the dashed curve 
must be repeated sequentially several times, and the majority of the measurement outcomes must be taken 
in order to correct errors (the repetitions are not shown). As in (a), logical or cTx correction operators 
(not shown) may be necessary on the output data block. 



Fig. |2|a); it consists of preparing a qubit in the state |+), measuring cr^ cr^ on this qubit and the input 
qubit, and finally measuring cf^ on the input qubit. 

The basic idea is now to encode the qubits in the repetition code. We have already mentioned that 
the preparation of a logical |+) involves simply preparing every qubit in the block in the state |+). A 
measurement of the logical cj^ is also very simple and robust; it can be implemented by measuring (j^ on 
each qubit in the block, and then taking the majority of the outcomes to correct errors. It remains to 
discuss how to implement a measurement of the logical on the two blocks. The difficulty is that the 
measurement must be robust against the dominant phase errors, and it must also respect the biased structure 
of the noise. Fig. [2|b) shows our solution: A single unencoded "ancilla" qubit is prepared in the state |+), 
and then CPHASE gates are applied between this ancilla qubit and all other qubits in the two blocks. Finally, 
a measurement of is performed on the ancilla, and a logical <Jx correction may be necessary on the output 
block conditioned on the measurement outcome. Since even a single phase error on the ancilla can cause an 
error in the measurement outcome, the measurement of the logical must be repeated sequentially 
several times; by taking the majority of the outcomes we can suppress the probability of a logical error at 
the output. 

Using the building blocks just discussed, it is possible to construct a fault-tolerant logical cnot gate that 
respects the physical-level bias of the noise; see Fig. [3j This circuit implements a logical cnot gate between 
the two blocks encoded in the repetition code, while at the same time the logical state of each block is 
teleported to a new block and phase errors are corrected. Because of the use of teleportation, this circuit has 
the additional feature that it largely prevents the propagation of leakage errors. In Fig. [3) the output data 
qubits always interact with an ancilla qubit prior to any interaction with the input data qubits; therefore, 
there is no possibility for leakage to propagate from the input to the output qubits. And furthermore, if we 
consider implementing a sequence of logical cnot gates, we observe that every qubit is eventually measured 
after only a small, fixed number of time steps, and the measurement effectively converts leakage to regular 
qubit errors [2T] . 

In addition, it is possible to construct with the same building blocks a fault-tolerant preparation of the 
logical |0) and |+) states, and a measurement of the logical <Jz and cf^ for the repetition code. The logical 
CNOT gate, the preparation of logical |0) and |+), the measurement of logical dz and ctx, and single-qubit 
measurements of exp(z^(Jz)crx for = 7r/4 and = tt/S suffice for implementing universal fault-tolerant 
quantum computation; the full scheme is discussed in our companion paper [19 . 

Since our building blocks only use CPHASE gates, |+) preparations, and measurements, and since phase 
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Figure 3: A logical cnot gate between two blocks encoded in the repetition code (here again, the 3-bit 
code). Ancilla qubits (shown in green) are prepared in the state |+), they interact via cphase gates with the 
data qubits (shown in red), and they are measured along the eigenbasis of a^. As in Fig.[2j the measurements 
with ancilla qubits must be repeated sequentially several times so that errors can be corrected by taking the 
majority of the outcomes (the repetitions are not shown). Finally, <Jx is measured on each qubit in the two 
input blocks, and the majority is taken on each block. Conditioned on the results of the majorities, logical 
correction operators (not shown) may be necessary on the output blocks |19] . 



errors commute with CPHASE gates, it is simple to estimate the error rates for the logical cnot gate given 
the physical-level error rates. In [19], we give upper bounds for the logical error rates in closed form as a 
function of the block size n of the repetition code and the number of repetitions k of the measurements 
executed with ancilla qubits. The outcome is qualitatively as in Eqs. ([T]) and ([2| with t = ck for some 
constant c ~ 2 or 3. It follows from this analysis that if the noise bias in our elementary operations is about 
10^ or greater, encoding in the repetition code is effective and logical errors are significantly weaker than the 
physical-level phase errors for s of order .1% or smaller; see Fig. [4][^ 

The noise bias reported in p!8l was about a factor of 10, which has motivated us to re-examine whether 
greater levels of bias are conceivable for the Koch qubit. As we will now discuss, the noise bias can, in fact, 
be dramatically improved at the expense of a very minor increase in the rate of phase errors. Fig. [5] shows 
the structure of the Koch qubit [14, 8 . It is nominally a "flux" qubit, meaning that the computational 
basis states |0) and |1) are quantum states corresponding to distinct circulating-current orientations; see 
Fig. |5j The mode of operation is highly tunable via an external "control flux" threading the small loop; 
see Fig. [6j At small <l>c, the degeneracy of the two circulating-current states is lifted by the "flux bias" 
e corresponding to the flux difference in the two large loops, and the two basis states are easily detected 
and initialized but not very phase-coherent. As <l>c is increased, the barrier decreases until it eventually 
disappears, and tunnel splitting between the basis states turns on rapidly. At even larger we enter an 
almost harmonic single- well regime, but before this another essential element of the Koch qubit enters the 
picture. 

The flux-qubit states are strongly coupled to the fundamental mode of a superconducting transmission- 
line resonator, so that |0) and |1) pass via an avoided level crossing to |0) and |i) corresponding to the 0- 
and 1-photon modes of the oscillator. The superconducting transmission line is highly phase-coherent [8 . In 
the IBM experiments we have seen circa 50,000 Ramsey fringes associated with these states corresponding 
to T2 = 2.5/isec [22]; we expect much longer T2 times to be possible. This fact has led us in [18] to the 
following strategy for implementing low-noise operations. When not being acted upon, quantum information 

^Ref. also discusses a more sophisticated procedure for decoding the repetition code which is effective for e of about .5% 
or smaher. For simphcity, we do not analyze this decoding procedure in this paper. 
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Figure 4: Upper bounds on the total probability °* of logical errors for the logical cnot gate as a function 
of the probability e of physical- level phase errors for a bias of 10^ and 10^. To obtain our bounds, we have 
optimized over the block size n of the repetition code, and the number k of repetitions of measurements with 
ancilla qubits; the optimal choice is to have n = while the optimal value depends on e and the bias as 
shown. The straight line with slope unity serves as a guide to the eye. 

is stored in the highly-coherent oscillator energy levels; we then say that the qubit is "parked." All needed 
operations are done by adiabatic pulsing out of parking. Each flux qubit has a fixed, untuneable two-qubit 
coupling to a set of nearby flux qubits in order to implement two-qubit gates. To assure that the effective 
coupling between parked qubits is negligible, these couplings are to be only between qubits with different 
resonator frequencies, so that resonant transfer of photons between different resonators does not take place 
in the parked state. Below we examine the simplest scheme with just two resonator frequencies, 3.1GHz and 
I • 3.1GHz (the commensurability of these frequencies aids in the maintenance of rotating frames for these 
qubits). Correspondingly, we have two species of qubits, A and B respectively, and two-qubit couplings exist 
only between qubits of different species — this is not a severe restriction since, in our encoding scheme, the 
only interactions are between data qubits (red qubits in the flgures) and ancilla qubits (green qubits). 

Fig.[7|^a) shows the energy levels of the coupled two-qubit system as a function of the control flux — assumed 
equal for the two qubits in the flgure — with the flux bias held on the "symmetric line" e = on which the 
effective qubit potential U has a reflection symmetry as in Fig. |6] In the "parking" region ^ lAG^o, the 
effective interaction between the two qubits is strongly suppressed. The energy levels are essentially those of 
the unperturbed transmission line resonators and, because to very high accuracy we have energy additivity 
^(|i)|i)) = ^(|0)|1)) +£^(|i)|0)), there is no conditional phase accumulation in the two-qubit state. Out of 
parking, this energy additivity is violated and a conditional phase shift can accumulate. If only one qubit is 
pulsed out of parking, single-qubit "diagonal" rotations of the form exp(— z^dz) can be effected (in the frame 
of reference rotating at the resonator frequency). If two coupled qubits are pulsed out of parking, a cphase 
gate can be implemented by choosing the pulse timing appropriately; at the same time, known single-qubit 
phase shifts accumulate which can be compensated by corrective diagonal rotations on each qubit [18 . 

There is, however, one point of concern which reveals itself more clearly when we plot the energy levels of 
the same two-qubit system as a function of time during the implementation of a cphase gate; see FiglT^b). 




At the energy- level crossing marked C, the avoided-crossing gap between the state and a state outside 

the computational space becomes small, and this may potentially lead to increased leakage errors. One 
possible strategy to avoid this problem would be to increase the avoided-crossing gap by changing the flux 
bias to e 7^ in order to move away from the symmetric line; this solution was adopted in [18 . However, 
departing from the symmetric line has its disadvantages. The lowering of symmetry makes the system less 
protected from low- frequency noise, and it also causes the effective relaxation time Ti to be shortened making 
it harder to achieve a high bias in the noise. 

It would therefore be desirable to perform the operation at the symmetric line. Our new observation is that 




■ 

Figure 5: Physical layout of two Koch qubits [8 coupled as needed to implement a cphase gate. Qubits 
A and B differ only in the resonant frequencies of the transmission line resonators attached to them. These 
frequencies are controlled by the physical length of the transmission lines, and are in the ratio 3:4. The single- 
qubit device Hamiltonians can be varied independently by controlling the "control flux" corresponding 
to the flux through the small loop, and the "flux bias" e corresponding to the flux difference in the two 
large loops. The SQUIDs perform quantum measurements on the states of the individual qubits. The qubits 
are coupled via a fixed mutual inductance M to a superconducting loop. For small values of the basis 
states |0) and |1) of each qubit correspond to different orientations, counterclockwise and clockwise, of the 
persistent current Ip. 

the avoided-crossing gap is actually very small; it is never larger than lOOkHz even taking parameter shifts 
due to low- frequency noise into account. So for our pulse profiles, the Landau- Zener tunneling probability 
is extremely close to one, meaning that we can pulse through the avoided-crossing gap and suffer essentially 
no leakage. After trying many pulse designs, and checking the noise bias they produce, we find the pair of 
pulses in Fig. [tJ^c) for the two species of qubits to work the best. We observe that in this pulsing scheme the 
low-frequency qubit B is unparked far more deeply (point Pb) than the high-frequency qubit A (point Pa); 
this choice is optimal because, if the two species of qubits were unparked symmetrically, the same level of 
dephasing noise could only be maintained by making the pulse duration significantly longer increasing the 
noise due to relaxation and leakage. 

To estimate the effect of noise in the implementation of our cphase gate, we have performed multiple 
computer simulations of the evolution of the two coupled qubits during the gate where, in each simulation, 
low-frequency noise is added to the ideal dynamics. To model low-frequency noise, in each simulation we 
take the actual applied flux and the actual timing in a pulse to deviate from the ideal by an amount which is 
chosen from a Gaussian distribution with zero mean and variance and 6psec respectively [18]. From each 
simulation we extract a unitary which describes the noisy implementation of the cphase gate; by averaging 
over a large number of simulations which corresponds to integrating over the Gaussian fluctuations, we 
obtain a superoperator describing the noisy implementation of the gate. By expressing this superoperator 
as A/" o CZ where CZ is the ideal superoperator when there is no noise, we can define an "error rate" in 
terms of the norm oi £ = M — X where X is a trace-decreasing superoperator proportional to the identity 
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Figure 6: The lowest energy levels of a single Koch qubit on the "symmetric line" corresponding to e = 0, 
as a function of ^c- For small <l>c = 1.4<^>o (*^o = h/2e)^ the ground state is doubly degenerate; the 
degeneracy may be lifted by setting e 7^ 0, allowing for qubit measurement and preparation. The two basis 
states |0) and |1) correspond to the two states \L) and \R) of a double-well potential with a very large 
potential barrier between them as sketched, which describe the two different orientations of the persistent 
current Ip in the flux qubit; the rescaled dynamical variable / is explained in [23 . For both |0) and |1), the 
superconducting transmission line is in its "vacuum" 0-photon state. As is increased to around 1.45<l>o, 
the barrier height drops, leading to a rapid increase in the tunnel splitting between the two lowest states. 
When the tunnel splitting equals the transmission-line resonator frequency cjt, there is an avoided level 
crossing with splitting 2g. For larger the qubit is "parked"; now, the two lowest energy states |0) and |i) 
correspond to the 0- and 1-photon states of the transmission line respectively, their energies are independent 
of <l>c, and their effective potential is highly harmonic as sketched. In this regime, for both |0) and the 
flux-qubit state is the symmetric state l^*) = + |^))- 

superoperator. Although this rate does not necessarily correspond to the probability of an error, it is possible 
to formulate a fault-tolerance analysis similar to the case of stochastic noise. In particular, we may expand 
^ = ^phase + Mother whcrc Sphase coutaius the tcrms of £ which are diagonal in the computational basis, and 
Mother contains all other terms. Then, Eqs. ([T]) and ([2| remain unchanged if e and e' are re-interpreted as 
the norms of ^phase and fother respectively. In the Supplementary Material, we give more details about this 
definition and about the derivation of the error rates that we discuss below. 

Based on our modeling, we find the following hierarchy of expected error levels for the cphase gate: The 
rate for phase errors is 1.96 x 10~^ for qubit A and 4.6 x 10~^ for qubit B. For all other types of errors, the 
dominant contribution is due to relaxation to the ground state during the 35.1nsec duration of the cphase 
gate; using our previous results based on the Caldeira-Leggett model which give an estimate of Ti ^ 10msec 
[23||T8], we expect their rate to be around 3.5 x 10~^ for both qubits A and B. Therefore, as desired, there is 
a large contrast between the rates for phase errors and all other types of errors. We must also finally assess 
the magnitude of leakage. Since leakage errors cannot in general be corrected by our repetition code, it is 
crucial that leakage is suppressed to a very large extent; in fact, we have found that the rate for leakage 
errors is approximately 3.5 x 10~^. 

We have also designed pulses for the preparation of a qubit in the state |+), for single-qubit diagonal 
rotations, and for the measurement of a^; these pulses are only slightly modified from [18]. Since the 
diagonal rotations commute with the cphase gate, they can always be moved to occur immediately before 
the measurements of a^] the combined operation is a measurement of an operator of the form exp(z^(Jz )crx. 
The rate for phase errors in a preparation of |+) is 2.75 x 10~^ for both qubits A and B, and the rate for all 
other types of errors is 3.5 x 10~^. However, because a preparation of |+) is performed by a non-adiabatic 
pulse, leakage is significantly larger than in the cphase gate; the leakage error rate is 3.77 x 10~^ for qubit 
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Figure 7: CPHASE gate modeling, (a) The energy levels of two coupled qubits for the special case = 

(b) Variations of the eigenlevels of the two-qubit system for the optimal implementation of a CPHASE gate. 
To very good approximation, the quantum evolution adiabatically follows the four states |0)|0), |0)|i), |i)|0) 
and |1)|1). A crucial moment in the evolution is at the crossing C (and its time-reversed image in the second 
half of the pulse), when crosses a state outside the computational space. Leakage to this state is 
suppressed because the avoided-crossing gap is very small (the inset shows the dispersion at C magnified 
by 200,000x), and because the angle of crossing is made larger by delaying the onset of the A pulse relative 
to the B pulse. Apart from this state and the four computational-basis states, all other eigenstates can be 
entirely left out of our numerical modeling because of their large distance from the and |i)|0) states. 

(c) The optimal control-fiux pulses for A and B. Pulses are constructed as sums of tanh functions in order 
to have a standard smooth shape. It is found preferable to unpark A much less deeply (to Pa) than B (to 
P5). Although Pa and P5 in (a) do not exactly correspond to the points Pa and P5 in (c) because and 
are not equal throughout the optimal pulses, they provide a good illustration of the relatively large energy 
difference (~ 1.5GHz) corresponding to the small change in (~ .004<l>o) between Pa and P5 in (c). 



A and 1.5 x 10~^ for qubit B. For a measurement of exp(i^(Jz)crx, we can describe noise in terms of effective 
errors of no specific structure preceding the ideal measurement, and their rate is 1.83 x 10~^ essentially 
independent of 0. 

Given these physical- level error rates, we can obtain upper bounds on the error rates for the logical cnot 
gate by using the equations in [19 . In this analysis, we use the error rates we have computed by averaging 
over a large number of simulations as fixed, constant error rates, and we will ignore any correlations between 
the fluctuations around these average values in space and time. The analysis must also consider leakage errors 
which are not discussed in [19 . Our method for analyzing leakage is based on the following observations. 
First, as we have already noted, the fact that teleportation is used to implement the logical cnot gate 
prevents leakage errors from propagating across multiple logical gates. In the worst case, one leakage error 
can propagate from the logical gate where it occurs to the following one. But this can be prevented by 
inserting a logical teleportation preceding every logical gate; i.e., the logical state of each output block from 
a logical gate is teleported to a new block as in Fig. |2] before the next logical gate is applied. Then, a leakage 
error that occurs in a logical gate or in the logical teleportations preceding it can only affect this logical gate 
and no other. 

In our analysis of the logical cnot gate, we have optimized over the block size n of the repetition code, 
and the number k of repetitions of measurements with ancilla qubits. For our physical- level error rates, the 
optimal choice is (n, /c) = (5,7), where k is larger than n because qubits of species B (which are chosen as 
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the ancilla qubits) are more noisy that qubits of species A (which are chosen as the data qubits). With this 
choice, we find that the logical cnot gate has nearly balanced rates for phase errors and all other types of 
errors of at most 4.62 x 10~^ and 3.98 x 10~^ respectively. 

This is an improvement by a factor of about 3 over the best alternative method we have for implementing 
a CNOT gate without the encoding in the repetition code; as we have already discussed, we have found no 
direct implementation of a cnot gate for the Koch qubit with an error rate better than 5%, and simulating 
the CNOT gate indirectly as in Fig. [l] leads to balanced rates for phase errors and all other types of errors of 
about 1.25%. Since our error-rate upper bounds also apply to the other logical operations for the repetition 
code needed for universality [19 , our analysis indicates that our estimated physical- level error rates for the 
Koch qubit are tantalizingly close to those needed for effective fault-tolerant operation; in the literature, 
the highest estimated error thresholds are of order 1%, and the best proven thresholds are of order .1% 

[ISlIiiliilEo]. 

To conclude, we have discussed an encoding scheme which has been especially tailored to the physical op- 
erations that are naturally available in the IBM qubit (cphase gates, qubit preparations and measurements), 
and the noise characteristics that can be expected for this system according to our theoretical modeling. 
We believe that the basic principles underlying our proposal — primarily, the implementation of adiabatic 
CPHASE gates with highly biased noise, and the suppression of phase errors by encoding in the repetition 
code — can be successfully adapted to apply to other promising systems besides flux qubits such as, e.g., 
superconducting phase qubits [H [7] . 

Our estimated error rates incorporate contributions from all sources of noise which are understood in 
experiments at present. However, the experimental reality today is that noise is dominated by Ti processes 
which are not fully understood, and coherence times are significantly below the values we have obtained 
from our calculation [18]. Our results should therefore be seen as preliminary and suggestive. Certainly, 
simplifications in modeling low-frequency 1/f noise have been made, so that noise fluctuations have been 
assumed to be constant during the execution of each gate, and also noise correlations and flux drifts across 
multiple gates have been ignored; see [18]. 

Furthermore, we should emphasize that there are several possibilities for complementing the methods of 
"software" error correction we have described here with "hardware" error correction which is done directly 
at the physical level; e.g., systematic noise correlations and flux drifts over long time scales could be sup- 
pressed by periodically recalibrating qubits off-line before they are re-used, or superconducting qubits could 
be designed to have physical error-correction properties as in [24 which uses the same set of elementary 
operations as ours. Finally, we should note that the topology of interactions for our repetition-code scheme 
is not attainable with short-range interactions on a square lattice. We expect however that a greater but 
limited range of interactions where each qubit is coupled to 10 or 20 other nearby qubits would suffice. Vari- 
ous possibilities exist for implementing such interactions for superconducting-qubit layouts where crossovers, 
multiple couplers, and indirect couplings via transmission lines are all available. 

DDV and BMT have been partly supported by lARPA under ARO Contract No. W911NF-04-C-0098. JP 
is supported in part by DoE under Grant No. DE-FG03-92-ER40701, NSF under Grant No. PHY-0456720, 
and NSA under ARO Contract No. W911NF-05-1-0294. 
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Supplementary Material on Error Rates 

p. Aliferis, F. Brito, D. P. DiVincenzo, J. Preskill, M. Steffen, and B. M. Terhal 



In our analysis, we model each noisy elementary operation by a superoperator. Therefore, we ignore tem- 
poral or spatial correlations between different operations, such as the correlations in time which are inevitably 
present for 1// noise. We consider two sources of noise. First, 1// noise in the control parameters — the 
applied fluxes, and the pulse synchronization — during the execution of each operation; 1// noise leads pri- 
marily to dephasing and it will be modeled by a superoperator A/i//. Secondly, thermal relaxation noise 
which is continuously present; relaxation is the primary source of errors different than phase errors and it 
will be modeled by an amplitude-damping superoperator A/ti • 

Leakage errors depend both on the implementation of our operations and also on 1// noise. We recall 
that during qubit preparation and measurement, which are implemented by non-adiabatic pulses, leakage 
can arise from Landau-Zener transitions to excited levels in the flux-qubit potential; during the execution 
of a CPHASE gate, leakage primarily occurs at the energy-level crossing marked C in Fig. ^h). Since the 
relevant parameter is the minimum energy gap between states inside and outside the computational space 
during the implementation of an operation, leakage errors can occur even in the absence of 1// noise. When 
1// noise is present, the minimum gaps are shifted from their ideal values when there is no noise, which in 
turn has an effect on leakage. 

We model a noisy preparation of a qubit in the state |+) as the ideal preparation followed by the noise 
superoperator 

AT^ATt, oATi/j , (3) 

where o denotes composition. We model a noisy CPHASE gate as the ideal gate followed by a superoperator 
as in Eq. (|3|, where A/i// is now supported on both qubit s acted upon by the gate, and A/ti is assumed to 
act independently on each of the two qubits. Finally, we model a noisy measurement of ex.p{iOaz)cr^ as the 
ideal measurement preceded by a superoperator as in Eq. ([3|. 

The noise superoperator for each elementary operation can be expressed in terms of a discrete set of 
Kraus operators (at most d'^ where d is the dimension of its support) [25 . If the identity is one of the Kraus 
operators, the noise model corresponds to local stochastic noise; in this case, we may define the phase error 
rate as the probability of all non-identity Kraus operators which are diagonal matrices in the computational 
basis. The probability of all other non-identity Kraus operators then defines the rate for all other types of 
errors. For this noise model, we can perform a fault-tolerance analysis to determine upper bounds on the 
probabilities of logical errors as in [19 . 

However, as we will discuss below, the superoperators in our modeling do not have this property so that 
noise cannot be simply described in terms of probabilistic errors. Nonetheless, we may use an alternative 
definition of an error rate, and with this definition the analysis in [19] remains essentially unchanged. The 
idea is to express JV as the sum of an ideal and an erroneous part, Af = X -\- where X is a trace-decreasing 
superoperator which is proportional to the identity superoperator [26]. We may then define a generalized 
error rate or error strength in terms of the distance between A/" and X, 

e^\\£\U = \\Af-i\U, (4) 

where || • ||o is the diamond norm [27 . If the superoperator £ has an n-qubit input, ||^||o = ||^n ^ ^||i = 
i^iaxx :||x||tr=i I K^n |tr, whcrc Xn is the identity superoperator on n qubits, and || • ||tr is the standard 
trace norm, i.e., ||A||tr = TrV A^A. 

In the remainder of this Supplementary Material, we give the details about the error rates for our el- 
ementary operations; the results are summarized in the end in Table [l] To simplify our calculations, we 
have estimated the various norms by only varying among a few possible inputs. Even though we have not 
performed a rigorous maximization, we believe that the inputs we have chosen are close to the worst case. 
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1// noise in cphase gates 

For the cphase gate, we have obtained the combined superoperator A/i// o CZ by integrating over the 
Gaussian fluctuations in the apphed fluxes and the pulse timing as described in the main text; here, CZ 
is the ideal cphase superoperator (up to diagonal rotations on each qubit which are all moved before the 
measurements and will be discussed separately). From these numerical simulations, we can then extract the 
Kraus operators for A/i//. Each Kraus operator is supported on a 16-dimensional space 
where Ti^ and Ti^ are 4-dimensional spaces corresponding to the two interacting Koch qubit s A and B. For 
both qubits A and B, 

7i«=7ig,,®Wga„s, (5) 

where H^^^ is a 2-dimensional space spanned by the flux-qubit states {\L)^\R)} and TY^^ns another 2- 
dimensional space spanned by the 0- and 1-photon states of the transmission line |^ Since prior to and 
after the implementation of the cphase gate information is stored in the transmission-line modes, the 
computational space corresponds to the tensor product of the two transmission-line spaces, H^fns — ^trans^ 
^5ans5 action of the gate on the space TY^^ = '^aux^^Aux ideally trivial, and any transfer of amplitude 
to this space corresponds to leakage. 

We have found that, within the precision of our numerical analysis only four Kraus operators carry 
significant weight while all the rest are negligible. These four Kraus operators {Mq, . . . , M3} are of the form 

Mfe = 4 + Mfc,rf + Mfe,^rf + Mfc,^ . (6) 

Here, Ik is proportional to the identity operator on 7Y^^, M^^d includes all terms that act as the identity 
on Hfl^ and as diagonal matrices in the computational basis on H^fns (gi^i^g ^ise to phase errors), M^^^d 
contains all remaining terms acting as the identity on TY^^ (giving rise to other types of errors in the 
computational basis), and finally Mk^i includes all terms that act non-trivially on TY^^ (giving rise to 
leakage errors). 

If we expand in the Pauh basis in H^^ , we find {/q = -9981 exp(a.2743)/^ (g) h = 0, /2 = 0, h = 0}, 
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(7) 



We may write A/i// = X + + E^d + ^d^ where i^(X) = XI^^q ^k^^l is trace-decreasing and proportional 
to the identity superoperator. Si contains all terms with at least one insertion of a leakage error {M^^i}^ i.e., 

3 

Si{X) = Miif{X) - Y,{h + + Mk^^d)X{h + + Mfc,.^)t , (8) 

S^d contains all remaining terms with at least one insertion of a non-dephasing error {M^^^d}^ 

3 

£^d{X) = {Miif - £i){X) - Y,{h + Mu4)X{h + Mfe,rf)t , (9) 

and £d contains all remaining terms with at least one insertion of a phase error {M/e ^^j, i.e., 

3 

£d{X) = {M^if -El- E^d){X) - hXll . (10) 



/c=0 



^ Since transitions to states with more than one photon in a transmission hne are neghgible, in our numerical simulations we 
truncate the infinite-dimensional space of the transmission-line modes at the first excited state; see |18| . 
^We use 8 decimal digits of accuracy. 
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We define the rate of leakage errors as the norm of £i] by taking as the worst-case input the state |1)|1) 
(c/. Fig.[7|;b)), we find \ \Ei\U ^ 3.5 x IQ-^. Similarly, we define the rate of non-dephasing errors as the norm 
of E^d] we find = O(10~^), which can be neglected since it is much smaller than the contribution 

due to relaxation to be discussed below. 

We finally define the rate for phase errors as the norm of Sd- By varying among several possible inputs, 
we obtained the largest value of ^ 4.73 x 10~^ for the Bell state 

|$o) = ^(|o)|o) + |i)|i)). (11) 

Here, |0) = |5')|0p) and |i) = \S)\lp), where \S) = + \R)) is the symmetric state in H^^^ and \np) is 

the n-photon state in Titans with Q either A or B; see Figs, [e] and [tI 

We may also estimate the rates for phase errors on qubit A and qubit B separately. For qubit A, we 

/f=T^^Ei^ S^d + £d^ 



modify our expansion by writing A/i// = -\- £i -\- S^d + where 



i^iX) = J2ih + M^^,)X{h + Mi^jt , (12) 



and M^^ includes those terms in M^^d which are proportional to (g) cj^ and so act trivially on qubit A. 
Then, £^ captures all terms that apply nontrivial phase noise to qubit A. By using the same Bell state 
as input, we find that phase errors on qubit A have a rate H^^f ||o ^ 1-96 x 10~^. If we perform a similar 
analysis for qubit B instead, we find H^^f ||o ^ 4.6 x 10~^. This shows that the effect of 1// noise on qubit B 
is stronger than on qubit A, which is physically expected since qubit A is unparked much less deeply than 
qubit B during the implementation of a CPHASE gate; see Fig. [7| 

1// noise in preparation 

For the preparation of the state |+), we have obtained the density matrix = A/i//(|+)(+|) by 

performing a similar integration over the fiuctuating fields as for the CPHASE gate. Depending on the species 
of qubit, /)|+^ is supported on li9 for Q either A or B, and it is of the form 

=^i+),d + ^i+),^ ; (13) 

here, is an operator supported on ^H^^ns ' ^1+),^ includes all remaining terms of Since 

after the preparation the information is stored in the transmission-line modes, the ideal state is |+)(+| = 
from which we can obtain by acting with phase errors alone; on the other hand, for 

all terms in /?|+),^ the state of the fiux qubit is orthogonal to l^*) so that leakage has occurred. 

We define the rate of leakage errors as £i = 1 17^|+)^/ 1 |tr; for qubit A we find £i ^ 3.77 x 10~^, while for 
qubit B we find ei ~ 1.5 x 10~^ (our pulses are not highly optimized to avoid leakage, and we believe 
the leakage error rate for qubit B can be improved if necessary). For phase errors, we define the rate as 
e = — c|+)(+|||tr where we are allowed to optimize over the choice of < c < 1; for both qubits A 

and B, we find e ^ 2.75 x 10"^. 

1/f noise in measurement 

A measurement of exp{iOaz)o-^ is implemented by applying the diagonal rotation ex.p{—iOaz) on H^ans^ 
followed by a non-adiabatic pulse mapping |+) and | — ) (where information is stored in the transmission-line 
modes) to |0) = \L)\Op) and |1) = \R)\Op) respectively (where information is stored in the fiux-qubit states), 
followed by a projection along the basis {\L)^\R)}. 

Diagonal rotations can be executed with pulses of short duration so that errors are very weak compared 
to other operations and can be neglected. We also assume that errors during the final projection can be 
neglected. The states \L) and \R) can be distinguished very accurately by setting very small (~ lA^o) 
so that there is a large potential barrier between them resulting in very high Ti in the {\L),\R)} basis {of. 
Fig. |6|; since the two states correspond to distinct circulating-current orientations, they induce different 
magnetic signals which can be detected by using the SQUIDs (c/. Fig.js]). 
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Qubit A {out = 27r X 3.1 GHz) 


Qubit B (c^T = 27r X 3/4 X 3.1 GHz) 


CPHASE 


£ 


i.yo X iu 






4.0 X lU 




e' 


3.5 X 10-^ 






3.5 X 10-^ 








3.5 X 10-^ 




1+) prep. 


£ 


2.75 X 10-^ 






2.75 X 10-^ 




£' 


3.5 X 10-'^ 






3.5 X 10-'^ 




£l 


3.77 X 10-'^ 






1.5 X 10-^ 


ex.p{i6az)(J^ meas. 


£ 


1.83 X 10-^ 






1.83 X 10-^ 



Table 1: Error-rate estimates for our elementary operations. For preparations and CPHASE gates, £ is the 
rate for phase errors, £' is the rate for all other types of errors, and £i is the rate for leakage errors. For 
measurements, the rate £ includes errors from all sources. 



To obtain the error rate for the remaining measurement process, we follow the evolution of the basis states 
1+) and | — ) by performing a numerical simulation similar to the cases already discussed. If the initial state 
is |+), we calculate the probability that the final state before the projection is orthogonal to |0), in which 
case we assume an error in the measurement outcome always occurs; and similarly for |— ), we calculate 
the probability that the final state is orthogonal to |1). We define this probability as the error rate for the 
measurement; for both qubits A and B, we find £ ^ 1.83 x 10~^. 

Relaxation 

We model relaxation noise by the amplitude-damping superoperator A/ti acting independently on each 
qubit; with Q either A or B, the two Kraus operators are 

Mo = + '-^^ , M, = ^a^il - af) , (14) 

where we have already defined and a^, and = I according to the tensor-product structure in 
Eq. ([5|. 

Since only Mi is non-diagonal in the computational basis, we define the rate for non-dephasing errors as 
1 1-^1 1 |o = 7, where Mi{X) = MiXM|. We also define the rate for phase errors due to the operator Mq as 
\\Mo-I\\<,, where Mo{X) = MqXM^ and i{X) = cX with < c < 1. If we take c = (l + vT^)^ /4, we 
find \\Mo-i\\o-l/2. 

For a CPHASE gate, we use the worst-case estimate Ti = 10msec, and we assume Ti can be treated as 
approximately constant during the execution of the gate; then, j = t/Ti = 3.5 x 10~^ where t = 35nsec is 
the duration of the gate. For diagonal rotations which are executed with pulses of duration t < 5nsec, 7 is 
very small and can be neglected. Finally, for preparation and measurement, Ti changes as a function of the 
control fiux and we calculate 7 = /q ^Stjt^ = 3.5x10"^. 
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